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1.  Introduction 


In  simple  terms  there  are  two  ways  of  fighting  wildland  fires:  (i)  directly 
attacking  the  combustion  with  water,  borate  bombs,  or  other  extinguishers,  (ii) 
indirectly  constructing  firebreaks  which  consists  essentially  in  removing  fuel 
along  certain  paths  of  some  width  which  circumscribe  the  fire  so  as  to  stop  further 
spread.  The  first  method  is  effective  for  early  suppression  of  the  development  of 
a  wildland  fire,  while  the  second  is  generally  required  for  well  developed  fires  of 
large  magnitude.  We  are  concerned  here  with  the  optimal  construction  of  firebreaks. 

Aside  from  the  greater  importance  of  constructing  firebreaks  since  they  apply 
to  larger  fires,  the  model  of  this  method  of  fire  fighting  lends  itself  more  easily 
to  an  evaluation  and  interpretation  of  the  parameters  involved,  particularly  when 
some  of  them  are  random  variables,  as  compared  with  a  model  for  direct  suppression. 
Further,  it  is  easier  to  describe  free  burning  fire  spread  models  for  the  study  of 
optimal  firebreaks. 


2.  Simplifying  Assumptions 

First  we  shall  assume  that  the  fuel  for  wildland  fires  is  distributed  uniformly 
on  a  flat  plane.  This  assumption  is  obviously  an  idealization  of  the  real  situation, 
but  it  may  be  possible  to  develop,  in  an  approximate  way,  projections  of  wildlands  on 
a  plane  so  that  equal  distances  in  the  direction  of  the  climatological  mean  wind 
vector  imply  equal  travel  times  for  a  free  burning  fire,  allowing  for  the  greater 
rate  of  travel  of  a  fire  up  slopes  as  against  the  velocity  of  spread  down  slopes  and 
also  for  variations  in  fuel  distribution.  Another  way  of  interpreting  this  assump¬ 
tion  in  the  real  situation  is  to  regard  the  idealized  velocity  of  fire  spread  in  the 
flat  plane  as  an  average  velocity  for  the  actual  wildlands.  In  any  event  one  cannot 
take  account  of  every  tree,  bush,  stream,  rock,  and  slope. 
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Next,  regarding  fire  spread  models  in  the  plane  for  free  burning  fires,  the 
fire  front  in  still  air  is  assumed  to  propagate  radially  with  constant  velocity  as 
illustrated  in  Figure  1.  If  a  mean  wind  vector  is  superimposed,  the  movement  of  the 
fire  front  is  assumed  to  be  modified  by  the  addition  of  a  fixed  velocity  vector  for 
the  center  of  the  fire  in  the  direction  of  the  mean  wind,  resulting  in  a  fire  front 
spread  as  illustrated  in  Figure  2.  Generally,  the  rate  of  propagation  in  still  air 
is  small  relative  to  the  wind  effect  and  the  contours  of  fire  spread  in  the  presence 
of  wind  will  be  long  and  narrow,  suggesting  an  approximation  in  the  form  of  a  plane 
wave  front  moving  with  constant  velocity  in  a  channel  of  fixed  width  L  as  illus¬ 
trated  in  Figure  3. 

The  details  of  analysis  presented  in  this  paper  will  be  restricted  to  the 
deterministic  plane  wave  fire  spread  model  of  Figure  3. 

A  study  of  optimal  strategies  for  firebreaks  related  to  the  fire  spread  models 
of  both  Figures  2  and  3,  with  wind  treated  as  a  random  variable,  will  be  reported  in 
a  doctoral  thesis  of  one  of  the  authors  of  this  paper,  to  be  published  later.  There 
a  periodic  review  policy  is  used  for  optimal  construction  of  firebreaks,  using 
dynamic  programming  to  determine  the  sequential  alteration  of  the  optimal  firebreak 
paths  for  arbitrary  sequence  of  realized  values  of  the  randomly  varying  wind. 

3.  Method  of  Constructing  Firebreaks 

Men  with  hand  tools  and  portable  equipment  are  used  to  construct  a  firebreak. 

A  total  crew  of  N  men  is  split  into  n  equal  groups  and  each  of  these  groups 
constructs  one  nth  of  the  firebreak.  For  n  even  or  odd,  the  firebreak  sections 
are  laid  out  symmetrically  in  a  continuous  path  as  illustrated  in  Figures  4  and  5. 
Starting  at  points  P  a  section  of  the  firebreak  is  constructed  to  a  point  Q  . 

Thus  the  fire  may  be  subdivided  into  n  contiguous  plane  wave  fires,  and  we  need 
only  study  For  given  total  crew  size  N  and  number  of  construction  groups  n  the 
optimal  firebreak  path  for  a  plane  wave  fire  of  width  ^  .  With  these  optimal  paths. 


Initial  Position 
of  the  Fire  Front 


0 


FIGURE  3.  PLANE  WAVE  FIRE  SPREAD 
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depending  upon  N  and  n  ,  one  may  determine  crew  size  and  group  size  so  as  to 
minimize  the  resulting  total  costs.  For  this  analysis  it  is  assumed  that  the  time 
required  to  distribute  the  fire  fighters  along  a  working  line  is  negligible. 

4.  Cost  Structure 

The  cost  structure  proposed  by  Parks  and  Jewell  [4]  is  used.  It  includes 
the  following  four  types  of  cost: 

(i)  A  fixed  cost  associated  with  maintaining  a  fire  fighting  organiza¬ 
tion  and  setting  it  into  action  at  the  time  of  a  fire,  denoted  by 

C  (in  dollars) . 

F 

(ii)  A  cost  proportional  to  the  crew  size  N  used,  including  items  such 
as  transportation,  and  other  "one-shot"  logistic  support  costs, 
denoted  by  C^M  where  Cg  is  dollars  per  man. 

(iii)  A  cost  proportional  to  the  total  number  of  man  hours  used  in 
constructing  firebreaks,  denoted  by  C  ’N*T  where  C  is  dollars 
per  man  hour  and  Tc  is  time  of  control. 

(iv)  A  cost  measuring  the  fire  damage  which  is  proportional  to  the  area 
burnt,  denoted  by  C  *A  where  C  is  dollars  per  acre  and  A  is 

D  D 

the  total  area  burnt. 

In  these  terms  the  total  cost  K  for  any  fire  attacked  by  a  crew  of  size  N 
is  given  by 

(1)  K  ■  CF  +  CS.N  +  Cm-N-Tc  +  Cb-A  . 

Although  not  explicity  indicated  in  (1)  ,  the  total  cost  depends  upon  the  path 

chosen  for  the  firebreak  and  the  number  of  construction  groups  n  ,  via  the  quantities 

T  and  A  .  In  order  to  see  this  we  need  to  consider  details  on  extremal  firebreak 
c 

paths  for  a  plane  wave  fire  of  width  ^  . 


FIGURE  4.  SECTIONALIZED  FIREBREAK  FOR  n 
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FIGURE  5 


SECTIONALIZED  FIREBREAK  FOR  n  - 
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5.  Extremal  Firebreak  Paths  (Given  N  And  n) 

Define  the  firebreak  path  by  a  real,  nonnegative  single  valued  function 
■  f(y)  ,  denoting  the  distance  of  the  firebreak  from  the  initial  position  of  the 
fire  front  for  all  values  of  y  in  the  interval  [0  ,  j^]  .  In  particular 
f(0)  =  Xq  >  0  (See  Figure  6). 

For  some  nominal  width  W  of  the  firebreak,  let  A  denote  the  firebreak  area 

constructed  per  man  per  unit  time.  Then  the  ratio  A/W  determines  the  velocity 

V  of  firebreak  construction  per  man  used, 
m 

The  position  of  the  fire  front  at  any  time  t  is  given  by 

(2)  xf(t)  *  Vf-t 

where  V^.  denotes  the  constant  velocity  of  movement  of  the  fire  front  and  the 
origin  of  time  is  chosen  so  that  t  =  0  corresponds  to  th*  Initial  position  of  the 
fire  front,  i.e.,  x^(0)  ■  0  . 

Denote  by  T(y)  the  time  at  which  the  fire  construction  group  reaches  the 
point  (f(y)  ,  y)  .  Then 

y 

(3)  T(y)  -  — 7JJ7  f  U  +  (f’(u))2]’s-du 

Vm'n)  0 

and  the  time  of  control  Tc  *  t(“)  is  given  by 

L/n 

T - 7J-  f  (1  +  (f'(u))2]’5-du  . 

m'n  /  0 

Thus,  the  cost  function  (1)  lepends  upon  the  firebreak  path  f(y)  and  the 
number  of  construction  groups  n  as  well  as  the  total  crew  size  N  .  Similarly, 
the  last  term  of  (1)  depends  upon  n  and  an  integral  of  f(y)  . 
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At  any  time  t  during  the  construction  of  the  firebreak,  the  fire  front 
cannot  pass  over  the  firebreak  path  f(y)  .  Hence  admissible  firebreak  paths  f(y) 
are  restricted  to  satisfy  the  constraints, 


(4) 


f(y)  > 


/ 


v  n . 

m\n  /  0 


[1  +  ( f  *  (u))2]^*du  V  y  c  [0  ,  -] 


In  these  variables  the  total  cost  function  (1)  becomes 


C  *N 
m 


L/n 

/ 


(5) 


K[f(y)  ;  n  ,  N]  =  C  +  C  -N  + 

V  - )  / 
m'n  /  0 


L/n 


[1  +  (f'(y))2]l5-dy 


l 


+  CB*n  |  f (y)dy 


and  extremal  firebreaks  are  defined  by  those  paths  f(y)  satisfying  (4)  which 
yield 

Min  K[f(y)  ;  N  ,  n]  ^ 

(6)  f(y)  V 

for  given  N  ,  n  . 


6.  The  Optimization  Procedure 

For  arbitrary  positive  crew  size  N  and  number  of  construction  groups  n  , 
extremal  firebreak  paths  are  determined  by  the  constrained  variational  problem  (6)  , 
in  which  the  endpoints  are  variable  on  the  boundaries  of  the  plane  wave  fire  front. 

In  this  minimization,  the  extremal  firebreak  paths  may  either  follow  the  fire 
front,  be  removed  from  the  fire  front,  or  be  a  mixture  of  these  two  possibilities. 

Having  determined  the  optimal  firebreak  path  as  a  function  of  N  and  n  , 
this  function  is  substituted  into  (5)  to  obtain  the  t  f  1  cost  as  a  function  of 
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N  and  n  ,  and  the  resulting  cost  function  is  then  minimized  with  respect  to  the 
number  of  construction  groups  and  total  crew  size. 

Thus  an  optimal  strategy  is  one  which  minimizes  cost  rel  tive  to  crew  size, 
number  of  groups,  and  firebreak  path. 

7.  Optimal  Firebreak  for  Given  N  and  n 

It  is  convenient  to  state  certain  general  propositions  concerning  the  extremal 
firebreak  paths  before  undertaking  the  variational  problem  (6)  . 

Proposition  1;  An  extremal  firebreak  path  may  follow  the  fire  front  if  and  only 
if 


v(-) 

m'n ' 


<  1  . 


Proposition  2:  If  — — -  <  1  and  the  initial  endpoint  of  the  extremal  firebreak 


m 


(!) 


path  satisfies  f(0)  =  xq  =  0  ,  then  the  extremal  path  entirely  follows  the  fire 

V, 


front.  If 


'(-) 

m'  n  f 


1  ,  there  is  no  feasible  control  with  this  endpoint  stipulation. 


Corollary:  An  extremal  firebreak  path  is  a  straight  line  follow-the-f ire 
front  path  starting  from  any  value  of  y  where  fire  front  and  firebreak  meet  on 
the  extremal  path. 

These  two  propositions  are  required  as  preliminaries,  because  the  standard 
variational  procedures  do  not  apply  for  a  straight  line  follow-the-fire  front  path 
which  permits  only  one  sided  variations. 

With  these  preliminaries,  there  are  evidently  two  parameter  situations  to  be 
considered  for  extremal  firebreak  paths, 


11 


(a)  p  =  7kp 1 

min  / 


(b)  0  =  1 
m\n/ 


Case  (a) :  p  >  1 


Here,  following  the  fire  front  is  not  admissible,  which  suggests  that  the 
constraints  (4)  may  be  simplified  to 


(7) 


L/n 

/ 


n  +  f  *  (y) 2  -  f '  (y)  -  0^- 1  dy  =  0 


In  fact,  define 


(8) 


G(y) 


f (y)  -  p yVl  +  f'(r)2dr 
0 


and 

(9)  G'(y)  =  f'(y)  -  p  +  f'(y)2  . 

If  p  >  1  ,  G'(y)  <  0  for  all  y  e  [0  ,  0  and  the  separation  between  men  and  fire 
front  is  a  strictly  decreasing  function  of  y  .  Thus,  condition  (7)  is  a  sufficient 
equality  replacing  (4)  ,  which  is  also  obviously  satisfied  when  f(y)  is  optimal. 
Hence  in  this  case  we  may  formulate  the  variational  problem  as  follows:  Define 


Yl 

(10)  J(f  (y))  *  f\Uy  ,  f(y)  ,  f'(y))  -  yj  -^yjf  1  dy 

yo 


f 
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where  yo  “  0  *  7i  =  ~  ,  Uq  is  a  multiplier  and 

(11)  L(y,f(y),f’(y))  =  A+f' (y)2  +  CBnf(y)  +  U0[p  ^l+f'(y)2  -  f  *  (y)  ] 

m 

Let  J(f(y))  have  an  extremum  for  f(y)  and  consider  variations  of  f(y)  defined 
by  f(y)  +  h(y)  ,  where  it  is  assumed  that  f(y)  and  h(y)  are  continuous  and 
differentiable  in  [0  ,  .  As  boundary  conditions  we  apply  6x^  =  h(yj)  , 

6xg  =  h(xg)  ,  allowing  the  endpoints  to  vary  on  the  lines  y  *  yo  *  y  =  yi  as 
illustrated  in  Figure  7.  For  arbitrary  h(y)  the  increment  in  J(f(y))  is  given 
by 

yi 

AJ(f(y))  *  J*  j  L(y ,  f  (y)  +  h(y),f'(y)  +  h'(y))  -  f  (y0)  +  6xq)  ]  j  dy 

yo 

yi 

-  f  jL(y,f(y),f'(y))  -Cy\ -y0')  f  (yp)  |dy 
yo 

and  the  corresponding  variation  SJ  is 


6J 


yi 

/ 

yo 


lLf(y) 


d 

dy 


Jf ’  (y)  I 


h(y)dy  +  L 


f  *  (y) 


6xi 


y*yi 


-  l 


f  ’  (y) 


+  u0 


6x0 


y=yo 


where  h(y)  ,  6xj  and  6xq 
f(y)  to  minimize  J ( f (y ) ) 

(12) 


are  arbitrary.  Hence,  as  necessary  conditions  for 
we  have 


Lf(y)  dy  Lf'(y) 


(13) 


Lf'(y) 


+  U0  =  0 


y=yo 


xq  +  <5xq 


x0 


FIGURE  7.  FIREBREAK  VARIATIONS  (p  >  1) 
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(14) 


Jf’(y) 


=  0 


y*yi 


Equation  (12)  is  the  standard  Euler  equation,  while  (13)  and  (14)  are 
transversality  conditions  for  the  variable  endpoints. 

Now,  using  the  definition  of  L(y  ,  f(y)  ,  f * (y) )  ,  we  obtain  from  (12) 


(12.1)  «  a(Un)  *y  +  B 

y l  +  f'(y) 

where  6  is  a  constant  of  integration  and 


(15) 


a(U0) 


C 

_2>  +  £ 

V  n 
m 


Uo 


Equation  (13)  becomes 


(13.1)  f'W  ?  =  0  =>  f •  (0)  ■=  0 

y  1  +  f’(0)Z 

and  0=0.  Equation  (14)  states 


(14.1) 


UO 


and  using  (12.1)  with  0=0  one  obtains 
(16)  U0  =  Cb*L 


and 


m  ,  r  L 

T  *  pCBn 
m 


(17) 


a(U0)  =  »o  = 
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Therefore  the  extremal  firebreak  path  for  case  (a)  is  given  by 


xm  =  f(y)  •  x0  I1  ' v/l  ~  (“oy)2i 


where  xq  =  f (0)  .  This  last  equation  may  be  written 


r  /  .  1  m2  .  2  1 

[x  -  (*o  +  — )  +  y  =  — r 

m  an  2 

a0 


to  show  that  the  optimal  firebreak  path  for  given  n  ,  N  is  a  circular  arc  with 
center  (xq  +  ~  ,  0)  and  radius  ~  ,  as  illustrated  in  Figure  8. 

The  initial  coordinate  xq  is  determined,  from  substitution  of  (18)  into 
(7)  ,  to  be 


-  iL.  sln-«+  i- 

an  \  n  /  an 


The  quantity  — —  is  less  than  unity,  since  p  >  1  for  case  (a)  and  aoy  <  1 

for  all  y  e  [0  ,  — ]  .  Thus  we  have  a  complete  solution  for  the  extremal  firebreak 

n 

path. 

One  final  remark:  the  cost  function  (5)  is  convex  in  f(y)  and  f'(y)  , 
and  the  necessary  conditions  (12)  ,  (13)  ,  (14)  are  sufficient  for  the  deter¬ 
mination  of  the  optimal  firebreak  path  for  any  given  positive  values  of  N  and  n  . 


Case  (b) :  p  <  1 


In  this  case  it  is  possible  for  the  firebreak  to  follow  the  fire  front 
starting  from  any  value  of  y  e  [0  ,  (see  Proposition  2  and  the  corollary 

following)  and,  in  order  to  investigate  the  composite  paths  which  may  occur,  let 


^  =  Min  y  3  f(y)  =  p 


1  +  f ' (r)  or 
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I 


Then  for  y  e  [0  , 
interval  [0  ,  y] 
limits  replaced  by 


(7.1) 


— ]  the  variational  problem  for  extremal  firebreaks  over  the 
is  one  which  minimizes  (5)  with  in  the  upper  integral 

y  .  As  in  case  (a)  we  shall  replace  the  constraints  (4)  by 


J  jp  A  +  f '  (y) 2  -  f '  (y)  +  j  = 


and  show  that  the  solution  satisfies  (4)  for  all  values  of  y  e  [0  ,  y]  .  The 

complete  path  of  the  optimal  firebreak  consists  of  this  solution  and  a  straight  line 

j  ^  dx 

with  initial  point  ( f (y )  ,  y)  and  slope  p/vi  -  p  ■  —  . 

Hence  for  case  (b)  we  use  the  definition  (11)  for  L(y  ,  f(y)  ,  f * (y ) )  and 
formulate  the  variational  problem  as  follows:  Define 


(21) 


J(f (y)) 


L(y  ,  f(y) 


f,<y>>  -  dy 


C  n  (y i  -  y) 

+  — T=T-  +  cb" 

m  V 1  -  p 


x(y i  -  y)  +  ^(yi  -  y) 


\/l-p2 


taking  a  firebreak  path  f(y)  for  y  e  [0  ,  y]  and  a  straight  line  follow-the-f ire- 
front  path  for  y  e  [y  ,  yj]  ,  where  yo  *  0  ,  yi  *  ~  ,  y  e  [0  ,  ^]  and  U0  is  a 
multiplier.  Let  J(f(y))  have  an  extremum  for  f(y)  (ye  [0  ,  y])  and  consider 
variations  of  f(y)  defined  by  f(y)  +  h(y)  ,  where  it  is  assumed  that  f(y)  and 
h(y)  are  continuous  and  differentiable  in  [0  ,  y)  •  The  boundary  conditions  are 


6x0  *  h(y0) 


6  k 


•<5y  =  h(y)  , 
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using  a  general  variation  (6x  ,  6y)  at  the  endpoint  (f(x)  ,  y)  as  illustrated 
in  Figure  9.  Then  for  arbitrary  h(y)  the  increment  in  J(f(y))  is  given  by 


y+6y 

AJ(f(y))  =  j  j  L(y,f  (y)+h(y)  ,f  '  (y)+h'  (y))  -  [f  (y0)-»-6x0  ]  |  dy 

yo 


J  j  L(y , f  (y) , f  ’  (y) )  -  f  (y0)  Jdy 


C  n 


+  i<yi“*~6,) '  (yi'y)S 

+  Cnn  j(*+6x)  (yi-y-6y)  +  ^(yi~y-6y) 2 -  *(yi-y)  -  ^(yi-y)2-/^ 


«••/  N/  1  J  ~  J  *  '  1  J  -  J  ’  f 

vl-p 

and  -he  corresponding  variation  6J  is 


6J 


/I 

yov 


Lf(y)  -  d7Lf’(y)  h(y)dy 


J  Lfl  ^  |  +  U0  (  6x0  +  , 

CBn(yi-y)  +  Lf,(y)  | 

(  y°yo  ) 

y-y 

L 


C  n 


y=y 


Xt 7' cBn(*+(y,‘^) 


-  Lr. 


f'<*> 


6y 


y-y 


where  h(y)  ,  6xq  ,  <5x  and  6y  are  arbitrary.  This  problem  is  one  of  Bolza  type 
[1]  .  Necessary  conditions  for  f(y)  to  minimize  J(f(y))  are: 


(22) 


L  —  — —  L  ■  0 

i (y)  dyLf'(y) 


,23) 


Lf • (y)  I  +  U0  "  0 


y=yo 
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(24) 


Lf » (y)  i  +  CBn(yi"?)  =  0 


y=y 


(25) 


C  n 

L  I  --S- 

y 


v  x/T”2" 

m  v  1-p 


-  CLnx  =  0 

D 


Here,  equations  (23)  and  (24)  are  transversality  conditions  for  the  variable 
endpoints  and  (25)  is  a  comer  condition  on  the  endpoint  with  general  variation. 

Now,  using  the  definition  of  L(y  ,  f(y)  ;  f'(y))  we  obtain,  as  before  in 
case  (a) , 


(22.1) 


\/l  + 


f'(y) 


f’(y) 


a(Uo)y  +  0 


Equation  (23)  has  the  same  form  as  (13.1)  and  implies  B  “  0  .  Equation  (24)  , 
although  modified  relative  to  (14)  ,  again  implies  U0  ■  C0L  and  a(Uo)  is  again 
defined  by  (17)  .  Hence  the  extremal  firebreak  path  for  the  interval  [0  ,  y]  is 
given  by  equation  (18)  . 

The  comer  condition  (25)  yields 


(26) 


£_ 


giving  a  determination  of  the  endpoint  (f(y)  ,  y) .  However,  since  p  <  1  , 
equation  (26)  does  not  surely  determine  a  value  for  y  which  is  less  than  . 

But,  using  (18)  for  the  extremal  path  in  the  cost  function  (5)  with  y 
replacing  ,  the  resulting  expression  is  convex  in  y  for  any  positive  N  and 

n  .  Hence  the  solution  for  y  may  be  written 


1  “0 

if 

anL  2  1 

p  <  or  n  <  | 

n 

R) 

V1 

if 

p  >  SLflii  or  > 

=  n  = 

R) 

(27) 


y 


1 
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The  initial  coordinate  x0  for  case  (b)  is  determined  (by  using  (18)  ,  (27) 
and  (7.1))  to  be 

(28)  x0  =  sin  1(p)  +  —  [  v/l  -  p2  -  1  ] 

a0  «0 

when  n  <  |a0L~^N  ,  otherwise  it  is  given  by  equation  (19)  . 

The  question  remains  -  does  this  solution  satisfy  the  constraints  (4)  . 

Using  (18)  in  (9)  ,  we  obtain 


*o(y  - 

G'(y)  =  7=^ . * 

VI  -  (a0y)Z 

and,  since  y  <  y  <  for  all  ye  [0  ,  50  ,  G(y)  is  strictly  decreasing  in  the 

—  UQ 

interval  [0  ;  ?]  to  a  zero  value  at  $  ,  Thus  the  constraints  (4)  are 
satisfied.  Moreover  aoy  <  p  <  1  for  all  y  e  [0  ,  y]  . 

The  general  form  of  the  optimal  firebreak  is  illustrated  in  Figure  10,  consisting 
of  a  circular  arc  connected  to  a  straight  line  follow-the-fire-f ront  piece  with  slope 
p/v/l  -  p^  ,  and  at  y  =  y  <  “  the  rlope  f'(y)  may  be  computed  from  (18)  and 
(27)  to  have  the  same  value.  Hence  the  condition  (25)  is  in  effect  a  Weierstrass- 
Erdman  corner  condition  [2]  . 


Summary  for  Cases  (a)  and  (b) 


In  Case  (a),  recall  that  p  >  1  and  <  1  . 


p  >  1  =>  n  > 


V  N 
m 

V, 


< 

n 


1  =>  n  >  agL 


Hence  it  follows  that 


FIGURE  10.  EXTREMAL  FIREBREAK  (p  <  1) 


ARBITRARY  CORNER  y 
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and 


0  V  N 

2  .  m 

n  =  aoL  — 


Thus,  if  we  regard  ?  to  be  a  variable  for  Case  (a)  as  well  as  Case  (b) ,  we  may 
state  for  both  cases  that 


(29) 


9  - 


0  V  N 

.  c  l  •  TT1 

if  n  <  a0L  — 


,  V  N 

j  r  »  f  m 

if  n  >  a0L  — 


and  the  optimal  firebreak  path  for  given  (N  ,  n)  is 


(30)  f (y)  = 


x0  +  ^  j  1  -  A  -  (a0y)2  { 


x„  +i-!i-x/TT72!+-?=£=f 


for  0  <  y  <  y 


for  V<y<— 
J  J  =  n 


where  p 


V 


m 


0 


and 


x0  =  £- 

«0 


-1  1  1  ’ 
sin  (a0y)  +  —  I  n  -  (a0y)‘ 

a0 


-  1 


(31) 


a0  = 


'B 


C.  ,  CBLVf 

V  V  N 
m  m 


8.  Optimal  Number  of  Construction  Groups  and  Optimal  Crew  Size 


For  arbitrarily  given  crew  size  N  ,  the  total  cost  under  optimal  firebreak 
path  depends  upon  the  number  of  construction  groups.  Treating  n  continuously. 
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we  investigate  the  cost  function  K(n  ;  N)  and  seek  the  optimal  number  of  construc¬ 
tion  groups  n*  for  arbitrary  N  >  0  . 

The  function  K(n  ;  N)  has  two  forms  depending  upon  whether 


Case  (i) 

Case  (ii) 


In  the  first  case,  it  follows  that. 


(32)  K(n  ;  N)  =  Cp  +  Cg- 


2  CR° 

+  - -  sin 


2(00)' 


i/oqL\ 

In/ 


and 


(33) 


_9K 

3n 


2(a0) 


sln-!(Mt)  .  (M) 


}>  0 


for  all  n  >  aoLl 


anL 


/V  N\ 

f  m  j 

\vf/’ 


since  this  inequality  implies 


°loL  < 


1  when  p  <  1  and 


<  1  for  p  >  1  (see  Case  (a)  above) .  Hence  K(n  ;  N)  is  a  monotone  increasing 

n  2  (wr 

function  of  n  for  all  N  in  the  range  n  >  aoLl 


t)- 


In  the  second  case,  it  follows  that 


(34) 


K(n  ;  N)  =  Cp  +  Cg 


V 


+  _ l 

\ v™  2n  7  7T 


2(oq) 


-1 


VI- 


=y  -  sin  (p) 


) 
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and 


(35) 


3K 

3n 


2  (  ot0  )  ' 


sin  1(p) 


-  ^  (M 

2C  v  1 

B  v  m  ' 


3/2 


(1  -  P  ) 


(36) 


92K  p 

CB  1  | 

fin) 

2  2 
(1  +  2p") 

3n2  n  A  -  p2 

_  2 (a0) 2  2CB  ' 

lVJ 

1  2  2 
(1  -  P  ) 

2  VmN 

for  all  0  <  n  <  agL— —  .  Although  (34)  is  not  monotone,  it  can  be  shown  that, 

f 

9^K  9K  9K 

since  p  <  1  in  this  case,  — r  <  0  whenever  —  =  0  and  —  *  0  .  Thus 

~  z  3n  9n 

3n  n=0+ 

/  2  VmN  \ 

over  the  range  I  0  <  n  <  agLr^—  1  the  function  K(n  ;  N)  is  quasi-concave. 

2  VmN 

Further  K(n  ;  N)  is  continuous  at  n  =  agL  — —  having  a  form  as  illustrated 

f 

in  Figure  11,  and  clearly  the  optimal  number  of  construction  groups  (not  restricted 
to  integers)  for  given  N  is 


(37) 


VI  V  N 


Then,  except  for  integer  disparity,  it  follows  from  (29)  that  the  optimal  fire¬ 
break  path  is  entirely  an  arc  of  a  circle. 

Turning  now  to  the  optimization  of  crew  size,  the  cost  function  K(N)  resulting 
from  (37)  is 


(38) 


K(N)  =  Cf+Cs.N+^ 


11  ~ 


a0LV{ 

V  N 
m 


V  N 

CB  V"°L  ~~v7 

- - -  sin 

2(a0r 


a0LVj 

V  N 
m 


FIGURE  11.  TOTAL  COST  AS  FUNCTION  OF  NUMBER  OF  CONSTRUCTION 
GROUPS  FOR  OPTIMAL  FIREBREAK  PATH  AND  ARBITRARY 


CREW  SIZE 
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where  otQ  depends  upon  N  as  given  by  (17)  .  It  is  convenient  to  make  a  change 
of  variable. 


(39) 


^  CgL 

Q  ~  V  N  =  C 

^  _ m  XT  i 


Then  the  cost  function  K(N)  is  written 


K(Q)  =  Cp  +  VfCBL(Q  -  1) 


m 


(40) 


it) 


7“  r3/2 

r  +  -7^ — — rr  sin 


'Q  -  1  (Q  -  1) 


W 


where  Q  is  a  monotone  increasing  function  of  N  and  1  <  Q  <  °°  as  0  <  N  <  °°  . 
Further 


(41)  K'(Q)  =  VfCBL 

m 


_  i  /js\  |  Q  +  1 
4  W  v/q(Q  -  1 


_+y/Q(3  -  (j)  -l/  1 

1)3/2  (Q.1)2Sin^ 


) 


Evidently  K (Q)  has  a  minimum  for  some  finite  value  Q°  >  1  ,  since  K(Q)  is 
positive  over  the  range  1  <  Q  <  00  with  K'(Q)  -»•  -<*>  for  Q  -►  1  and  K'(Q)  >  0 
for  Q  -*■  00  . 

Hence,  treating  N  continuously  in  the  range  0  <  N  <  00  ,  the  minimum  of 
K(N)  is  obtained  at  a  value 


(42) 


N°  =  (Q°  -  1)CrL 


(c  ) 

\  m  / 


and,  if  the  maximum  available  number  of  men  is  N  the  optimal  crew  size  is  given  by 

l  N°  if  N°  <  N 


N*  = 


u 


N  if  N°  >  N 


(43) 
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With  this  optimal  value  of  N  ,  the  optimal  number  of  construction  groups  n*  is 
determined  from  (37)  and  in  turn  the  optimal  firebreak  path  is  computed  from 
equations  (29)  ,  (30)  ,  and  (31)  • 

Note  that  this  solution  always  yields  a  circular  arc  firebreak  meeting  the 
fire  front  during  construction  only  at  the  endpoint  y  =  “  .  But,  if  we  restrict 
N  and  n  to  be  integers  such  that  j  is  an  integer,  the  optimal  firebreak  may 
have  a  terminal  straight  line  portion  following  the  fire  front. 

In  order  to  obtain  the  integral  optimal  values  for  crew  size  and  number  of 
construction  groups,  define 

[a]  =  largest  integer  equal  to  or  less  than  a 
<a>=  smallest  integer  equal  to  or  greater  than  a  . 


Then  compare  the  following  four  costs: 


(i)  K  [n*]  ,  [n*] 


N* 

n* 


(ii)  K  ^[n*]  ,  [n*] 


(iii)  K  (<n*>  ,  <n*> 


N* 

n* 


(iv)  K  (<n*>  ,  <n*> 


and  select  the  corresponding  integral  optimal  values  for  n  and  N  . 

Notice  that  the  following  subproblems  which  may  be  of  interest  in  themselves 
have  also  been  solved: 

(1)  Given  N  and  n  ,  finding  an  optimal  firebreak  path. 

(2)  Given  N  ,  finding  the  optimal  number  of  construction  groups  and  the 
optimal  firebreak  path  for  each. 

(3)  Given  that  N  <  N  ,  finding  the  optimal  values  of  N  and  n  as 
well  as  the  optimal  path  for  the  firebreaks. 


I 
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It  is  interesting  to  observe  that  certain  intuitive  results  for  particular 
cases  are  obtained  by  substituting  the  relevant  values  of  the  parameters. 

For  example,  when  Cg  =  0  (i.e.,  there  are  no  transportation  or  "one-shot" 

logistic  costs)  and  unlimited  men  are  available,  we  obtain  N*  =  00  ,  n*  =  «>  , 
xq  =  0  and  the  optimal  firebreak  consists  of  a  straight  line  which  coincides  with 
the  initial  position  of  the  fire,  recalling  that  we  assumed  the  time  required  to 
get  the  'ire  fighters  distributed  along  a  working  line  was  negligible.  Moreover 
the  related  optimized  total  cost  is  finite,  being 

C 

C  +  L  — — 
lf  +  L  V 

m 


obtained  by  computing  lim  K(Q)  from  (40)  . 

Q-x» 

Next  consider  =  0  ,  i.e.,  a  situation  where  fire  fighting  is  done  by 
volunteers.  Then,  it  follows  from  equations  (33)  and  (35)  that  n*  -  1  for  any 
N  ,  and  from  (40)  it  follows  that  N*  ■  1  .  Then  from  (29) 


9 


if 

if 


Vc  <  V 
f  m 


V  >  V 
f  =  m 


and  optimal  firebreak  paths  of  both  types  arise  depending  upon  the  relative  size 

of  V£  and  V  . 
f  m 

If  C_  =  0  ,  it  is  clear  from  the  structure  of  the  problem  that  N*  =  0 

D 

implying  that  it  is  optimal  not  to  fight  the  fire. 

9.  Stochastic  Models  and  Related  Optimization 


So  far  we  have  considered  only  a  deterministic  model  for  the  study  of  an 
optimal  strategy.  The  most  significant  random  variable  to  be  allowed  for  is  the 
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velocity  of  fire  spread  .  Taking  as  a  nonnegative  random  variable  with 

known  probability  distribution,  it  is  reasonable  on  practical  grounds  to  assume  a 
finite  upper  bound  on  the  realized  values  of  this  random  variable.  Then  an  optimal 
strategy  is  one  which  minimizes  expected  cost. 

For  this  extension  a  multi-stage  dynamic  programming  decision  model  may  be 
used,  in  which  the  dispatching  of  firebreak  construction  units  is  done  periodically 
based  on  observed  values  of  the  velocity  of  fire  spread  at  predetermined  points 
of  time. 

As  in  the  deterministic  model,  a  number  of  construction  groups  n  is  used 
with  an  assumption  that  each  one  nth  of  the  fire  front  has  the  same  realized 
value  of  fire  velocity  at  the  beginning  of  any  review  period.  The  value  of  n 

is  taken  fixed  for  all  periods  in  order  to  determine  the  optimal  policy  for  any 
number  of  construction  groups,  and  the  resulting  suboptimized  cost  function  K(n) 
may  be  minimized  to  determine  the  optimal  number  of  construction  groups  n*  . 

For  arbitrary  n  ,  the  optimal  policy  is  a  feedback  control  rule  which 
minimizes  the  total  expected  cost  for  an  unknown  number  of  periods,  in  which  an 
equal  number  of  men  are  dispatched  to  each  construction  group  at  the  beginning  of 
each  review  period  based  upon  the  realized  values  of  fire  velocity  in  all  previous 
periods.  In  practice  the  optimal  firebreak  path  during  any  period  is  determined 
in  accordance  with  the  realized  value  of  fire  velocity  at  the  beginning  of  the 
period,  and  with  the  terminal  point  of  the  firebreak  for  the  previous  period  as  the 
initial  point  for  the  firebreak  in  the  current  period.  The  policy  structure  is 
open  end,  i.e.,  no  fixed  number  of  review  period  are  set  in  advance,  and  the 
terminating  stage  is  reached  when  the  men  enclose  the  fire  front.  For  the  extremal 
problem,  the  main  difference  relative  to  the  deterministic  case  is  that  the  vari¬ 
ational  problem  has  a  fixed  initial  point  for  all  periods  except  the  first.  The 
initial  point  of  the  firebreak  in  the  first  period  is  a  free  endpoint  of  the 
related  variational  problem  and  it  is  determined  by  the  optimization  of  cost  with 


31 


respect  to  n  . 

If  it  arises  that  the  fire  front  overtakes  the  fire  line,  then  the  policy 
operates  to  dispatch  additional  men  to  each  construction  group  so  as  to  construct 
a  follow-the-f ire  front  firebreak.  This  involves  a  cost  of  dispatching  in  addition 
to  the  operating  cost  of  the  extra  men,  and  the  optimal  policy  is  determined  under 
this  rule  of  dispatching.  Since  the  transportation  cost  Cg  is  incurred  the 
instant  a  man  is  sent  and  the  total  cost  is  otherwise  a  decreasing  function  of  the 
construction  group  size,  it  is  never  optimal  to  call  men  back  until  completion  of 
the  firebreak. 
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